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We consider the effects of polydispersity on isotropic-nematic phase equilibria in thermotropic 
liquid crystals, using a Maier-Saupe theory with factorized interactions. A sufficient spread (« 50%) 
in the interaction strengths of the particles leads to phase separation into two or more nematic 
phases, which can in addition coexist with an isotropic phase. The isotropic-nematic coexistence 
region widens dramatically as polydispersity is increased, leading to re-entrant isotropic-nematic 
phase separation in some regions of the phase diagram. We show that similar phenomena will occur 
also for non-factorized interactions as long as the interaction strength between any two particle 
species is lower than the mean of the intra-species interactions. 



I. INTRODUCTION 

Solutions of anisotropic colloids or stiff polymers can 
exhibit phase transitions from isotropic to orientationally 
ordered nematic states as temperature is lowered or den- 
sity is increased. Such systems are often polydisperse, 
containing e.g. colloid particles with an effectively con- 
tinuous range of lengths and diameters or polymers with 
a distribution of chain lengths. Such polydispersity has 
experimentally measurable consequences. For example, 
in polydisperse colloidal rods with nearly hard interac- 
tions, demixing into two nematic (N) phases and coex- 
istence of two nematics and an isotropic (I) phase can 
occur P, Q • This is to be contrasted with the hypotheti- 
cal case where all colloidal particles are identical, i.e. the 
solution is monodisperse; there at most two phases (I— N) 
can coexist. 

For the theoretical description of phase transitions to 
nematically ordered states there are two broad classes of 
models. The first focusses on short-range repulsions be- 
tween particles and is therefore appropriate for lyotropic 
liquid crystals whose phase changes are controlled pri- 
marily by density, ft contains in particular the Onsager 
model of hard rods, whose second- virial truncation be- 
comes exact in the limit of long thin rods Q, and the 
Flory lattice model Q. The second class of models 
concentrates on the anisotropy of longer-range attrac- 
tive interactions, ft contains in particular the Maier- 
Saupe model, a mean-field theory which expands the 
orientational interaction dependence in low-order Legen- 
dre polynomials 0, @- Models of this type are useful 
for thermotropics, where temperature is the dominant 
control parameter. Theories combining these two ap- 
proaches have also been developed, see the review Q, 
and extended to include e.g. particle flexibility @. 

For lyotropics, binary (and to lesser extent ternary) 
mixtures have been thoroughly investigated theoretic ally 
using the Onsager and Flory models jj| 0, 0, 0, 
IT^Il^lTrilT7lfl8( . and phase separation into two nematic 
(N-N) phases and three-phase I-N-N coexistence pre- 
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dieted for sufficiently disparate rod sizes. Early work on 
binary thermotropic mixtures ignored the possibility of 
fractionation, i.e. unequal partitioning of different parti- 
cle species across coexisting phases. Within this approx- 
imation a mixture behaves like a monodisperse system, 
with no I-N coexistence region 0, |5(J. Later studies 
that accounted for fractionation effects did find I-N |2l| 
and N-N j2^| coexistence in bidisperse Maier-Saupe the- 
ory. We note as an aside that similar effects have also 
been seen in theories of mixtures of semi-flexible or stiff 
polymers; see e.g. H, H3, 0, E3 ■ However, there a driv- 
ing force for demixing even between isotropic phases is 
normally present, e.g. because of unfavourable interac- 
tions between unlike species. N-N demixing then occurs 
rather trivially when the nematic order prevalent at low 
temperatures is 'superimposed' on this isotropic demix- 
ing transition. 

Theoretical studies of genuinely polydisperse systems 
have been rather rarer. For distributions of the polydis- 
perse attribute (particle length, diameter etc) with two 
narrow peaks one of course expects phase behaviour sim- 
ilar to the bidisperse case [2(J H3, but it is not at all 
obvious what happens for unimodal (single-peaked) dis- 
tributions. For lyotropics, work on the polydisperse Fl ory 
model [2^, |2ii 133, El focussed on general features [32j 
such as fractionation and the widening of the I-N coex- 
istence region, but did not investigate the possibility of 
N-N demixing. For the Onsager model, primarily poly- 
dispersity in particle lengths has been studied. For a 
narrow distribution of lengths, perturbative calculations 
show that the I-N coexistence region again widens [33| . 
Recent studies [3 HH f° r wider length distributions have 
confirmed this and provided evidence for a region of I 
N-N coexistence in the phase diagram. However, the 
existence of N-N demixing at higher density remains an 
open question; in simplified Onsager theories N-N coex- 
istence is found only in a limited density range |36l l3?| 
or not at all pf . 

For polydisperse thermotropics, Sluckin used Maier- 
Saupe theory to predict the opening up of an I-N coex- 
istence region in slightly polydisperse systems |39] . This 
was based on a perturbation theory around the monodis- 
perse case, which did not allow for a study of possible N- 
N demixing. Semenov obtained similar results for ther- 
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inotropic semi-flexible polymers with small polydisper- 
sity in the chain lengths [4(j . 

Our aim in this paper is to find out whether a unimodal 
distribution of the polydisperse attribute is enough to 
cause N-N demixing and the associated I-N-N coexis- 
tence in polydisperse thermotropics. We use the simplest 
Maier-Saupe theory for this purpose, with factorized in- 
teractions. As explained in the next section, this gives a 
truncatable free energy for which phase equilibria can be 
calculated relatively efficiently. The resulting phase dia- 
grams are presented in Sec. IIIII and wider implications 
and possible generalizations are discussed in Sec. IIVI 



II. METHOD 

We consider a system with liquid crystal particles hav- 
ing a continuous range of values of some polydisperse 
attribute I. This could be particle length, for example, 
or polarizability in the case of van der Waals interactions; 
for definiteness we refer to length below. The state of a 
single phase of a system is then described by a distri- 
bution n(l, f2), defined so that n(l, ft) dl dVt/A-K gives the 
fraction of particles with length in an interval dl around /, 
and orientations O in a solid angle rffi. This distribution 
is the natural extension for a polydisperse system of the 
usual orientational distribution P(f2) used for monodis- 
perse systems of rods ^lj . The rod orientation Vl can be 
parameterized in terms of the angle 9 with the nematic 
axis and an azimuthal angle ip; due to the cylindrical sym- 
metry of the nematic phase n(l, fl) is independent of ip. 
Using c?Sl = dcosO dip the length distribution, obtained 
by integrating over orientations, is therefore 



dtt p(l,n) = / d6n(l,0) 



(1) 



where Jdd . . . = (1/2) J , dcosO .... The orientational 
distribution of the particles can be factored out from 
n{l,9) as H3 



n(l,0) =n{l)P l (6) 



(2) 



where Pi (9) represents the probability of finding a rod of 
given length / in orientation f2 = (9, ip) and is normalized 

according to Jd9Pi(9) = 1. In the isotropic phase, one 
has Pi(8) = 1 and n(l,9) = n{l). 

In the above notation the free energy per particle of 
polydisperse Maier-Saupe theory [39^ can be written as 
/ = /id + /ox, where the ideal part /id is the free energy 
of an ideal mixture while the excess part is 

f e * = ~fdl dl' n(l)n(l')u(l, l')S(l)S{l') (3) 

Here S(l) = Jd9 Pi (9)V2 (cos 9) is the orientational or- 
der parameter of particles with length I, defined as the 
average of the second Legendre polynomial 7-2 (cos 6?) = 
^(3cos 2 # — 1). The coefficients u(l,l') determine the 



strength of the attractive interaction favouring nematic 
ordering and depend on the particle lengths / and I 1 . We 
note that the description of the system in terms of a nor- 
malized distribution over lengths and orientations rather 
than a density distribution contains the implicit assump- 
tion that density variations with temperature or between 
coexisting phases are weak and can be neglected. Ac- 
cordingly, we will not concern ourselves with the density 
dependence of the interaction strengths u(l,l'). Maier 
and Saupe [j| assumed this to be quadratic, but Cot- 
ter ^3 later showed that a thermodynamically consis- 
tent derivation of Maier-Saupe theory as a mean-field 
approximation requires a linear density scaling. 

In the following we make the common assumption that 
the interaction strengths factorize according to u(l, I') = 
[u(l, l)u(l', V)} 1 / 2 . It is then sensible to switch from I to 
u ee u(l, I) as the polydisperse attribute, giving for the 
ideal and excess free energies (fee = 1) 



/id = T dud9 n(u,9)\nn(u,9) 



T / dun(u) 



1 



lnn(w) + / d9P u (9) In P u (9) 



(4) 



(5) 



/ ox = — — / dudv! n(u)n(u')Vuu' S(u)S(u') (6) 

With our factorized interactions the underlying nature of 
the polydisperse attribute has become irrelevant; poly- 
dispersity manifests itself only through the spread in the 
interaction strengths u. 

For a given distribution n(u), the orientational dis- 
tributions are obtained by minimization of the free en- 
ergy with respect to the orientational distributions P u {9). 
Bearing in mind that these are normalized, one finds 



Pu(0) 

z{u) 



V) exp (/^v^) 
d9 exp (/3Sy/uP 2 ) 



(7) 

(8) 



where we have abbreviated j3 = 1/T and V2 = ^(cosfl) 
and z(u) is the normalizing partition function for P u (9). 
The order parameter that appears in I|7I8II is 



dun(u)\fuS{u) — dud9 n(u,9)^/uV2 (9) 



s 



and determines the excess free energy / cx = — S 2 /2. It 
obeys the self-consistency equation 

S= fduTOniu)^^^ (10) 



/ d6 exp (j3S^V 2 ) 



This always has the trivial solution S = corresponding 
to an isotropic phase; nematic phases are characterized 
by S > and become physically relevant where they have 
a lower free energy than the isotropic solution. 

To calculate phase equilibria we need the chemical 
potentials n(u), which need to be equal in coexisting 
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phases. Denote N(u) the particle number distribution, 
so that the total particle number is N — J duN(u), and 
n(u) = N(u)/N. Then = SF/SN(u), where the 

extensive free energy is F = N f[N(u)/N]. This gives 



Sf 
Sn(u) 



+ /- / du'n{u') 



Sf 



8n{u') 



(11) 



In evaluating the derivatives Sf/Sn(u) we do not need to 
take into account the variation of the P u {9) because they 
are chosen to minimize /. A short calculation then gives 
6f/6n(u) = Tkx(n(u)/z(u))+T and so 



fj,(u) = Tin 



n(tt) 
z(u) 



(12) 



If the system separates into P phases a = 1, . . . , P, their 
it-distributions n a (u) can therefore be written as 



n„(u) = R(u)z a (u) exp(-/3^/2) 



(13) 



The common factor R(u) = exp[/?/Lt(u)] follows from the 
requirement of particle conservation: if phase a contains 
a fraction x a of all particles, then ^2 a x a n a {u) = no(u), 
where no (it) is the overall or "parent" distribution of 
u. Using also (JJJ, the joint distribution over interaction 
strengths and orientations in each phase therefore takes 
the form 



i a {u,9) 



"o 



(u)exp(/95„VuP2-/9^/2) 



J2b x b J dO exp (flSby/uPz 



PSI/2) 



(14) 



1 = / dudd n a (u,9) 
S a = dud9 n a (u,9)y/uV2 



A phase split involving P phases is therefore character- 
ized by the 2P variables S a and x a ■ These are determined 
by as many conditions, namely 



(15) 
(16) 



The relations (|15fl automatically guarantee that ^ a x a = 
1 because the parent no (it) is normalized by assumption. 

In our numerical work, we exploit the fact that the ex- 
cess free energy f cx — —S 2 /2 is truncatable 0, i.e. it in- 
volves only a single moment S of the distribution n(u, 9). 
Phase equilibria can then efficiently be found using the 
moment free energy (MFE) method |||4||4|. This 
constructs a MFE which depends only on S, and from 
which phase coexistence can be predicted in the standard 
manner. In this simplest form the MFE gives only the 
onset of phase separation exactly. However, by including 
additional moments, defined by adaptively chosen weight 
functions, increasingly accurate solutions are obtained. 
Using these as initial points, we are then able to find the, 
for our free ene rgy, exact solutions of the phase equilib- 
rium equations [361 l3?l l3c| . We omit the details of the 
implementation because they are similar to our study of 
an approximate Onsager model |36ll37j| . Differences arise 
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FIG. 1: Phase diagram for systems with a Schulz distribu- 
tion of interaction strengths u. The nature of the coexisting 
isotropic (I) and nematic (N) phases is shown as a function 
of the inverse temperature and the parent polydispersity S. 
Dashed lines indicate the expected continuation of the phase 
boundaries where they cannot be determined reliably from 
the numerical calculations. The inset shows the whole range 
of polydispersities studied, including the monodisperse limit 



primarily because, in the Maier-Saupe model, we are de- 
scribing the system in terms of normalized distributions, 
as would be appropriate in other systems for the case of 
fixed pressure rather than fixed volume; the appropriate 
modifications of the MFE approach are described in • 



III. RESULTS 

We concentrate on the example case of a Schulz parent 
distribution of interaction strengths, given by 



no(it) oc u z e 



-(z+l)u 



(17) 



The coefficient in the exponent has been chosen such that 
the mean interaction strength (u)o is fixed to unity, thus 
setting our energy and temperature scale. The width of 
the distribution is characterized by its standard deviation 
divided by its mean. This quantity, denoted 5 below, is 
often also referred to simply as the polydispersity and 
is related to z by 6 = (1 + z)^ 1 / 2 . To avoid interac- 
tion strengths which are arbitrarily large compared to 
the mean, we impose an upper cutoff of u max = 4. This 
affects the mean value of u only negligibly in the range 
of 5 of interest to us, while the relation between S and z 
has to be calculated numerically. 

Fig. n shows the resulting phase diagram, indicat- 
ing for each given parent polydispersity 5 (in the range 
S = . . . 0.60) and temperature the nature of the phases 
coexisting at equilibrium. Focussing first on the inset, 
which shows the whole range of 6, we see that in the 
monodisperse limit 5 — > there is no coexistence gap, 
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FIG. 2: Phase diagram for a bidisperse mixture with inter- 
action strengths ui = 1 and U2 = 3.5; x is the fraction of 
particles of the second kind in the parent. The dashed line 
indicates the temperature where three phases (I-N-N) can 
coexist. 



FIG. 3: Distributions of particles according to interaction 
strength, n(u). The solid lines show the three coexisting 
phases (I-N-N) at 5 = 0.571 and 1/T = 15; the dashed line is 
the parent distribution. Inset: Same for the I-N coexistence 
at S = 0.571 and 1/T = 5.5. 



with the transition from the isotropic (I) to the nematic 
(N) phase taking place at the well-known Maier-Saupe 
temperature of 1/T C ps 4.54. At nonzero S an I-N coexis- 
tence region appears. Its width initially grows quadrati- 
cally with 5, as predicted by perturbation theory for small 

5 39]. However, as 8 increases, the two boundaries of 
this I-N region are affected in an increasingly asymmet- 
ric way. The onset temperature of nematic ordering com- 
ing from high T changes rather little, increasing by only 
around 12% as the polydispersity increases from zero to 

6 = 0.6. The lower boundary of the I-N region, on the 
other hand, is shifted to significantly lower temperatures. 
In fact the broadening is so strong that the phase bound- 
ary becomes re-entrant: for 0.553 < 5 < 0.567 one has 
the phase split sequence I — > I-N — > N — > I-N — > ... as 
temperature is lowered. 

For polydispersities above S » 0.45, Fig.[I]shows in ad- 
dition that the system can fractionate from a single-phase 
nematic into two and later even three coexisting nemat- 
ics as T is lowered. The boundary of the N-N region 
meets that of the I-N region, and at that point a three- 
phase coexistence region (I-N-N) opens up where one 
isotropic and two fractionated nematic phases coexist. 
The structure of the phase diagram in this area is similar 
to the high-density phase diagram of polydisperse hard 
spheres l47_ . where two or more fractionated solids 
can occur with our without a coexisting fluid phase. 

It is worth emphasizing that all coexistence regions 
in the phase diagram are caused by polydispersity; the 
monodisperse system only has an I-N transition at a sin- 
gle temperature. Also, the presence of an effectively infi- 
nite number of different particle species, characterized by 
their respective interaction strengths u, means that the 
number of possible phases can be in principle arbitrary. 



This should be contrasted with the case of a binary mix- 
ture _2_J , which is illustrated in Fig. [__ for particles with 
interaction strengths u\ = 1 and u-i = 3.5; x is the frac- 
tion of particles of the second kind in the parent. A 
region of N-N demixing again occurs at low tempera- 
tures, but I-N-N coexistence is possible only at a single 
temperature. 

The rich phase behaviour that occurs in the polydis- 
perse system is driven by fractionation, i.e. by the fact 
that particles of different interaction strength u partition 
themselves unevenly among coexisting phases. As an ex- 
ample, we show in Fig. [31 the distributions n a (u) for the 
three coexisting phases (I-N-N) into which a Schulz par- 
ent with 5 = 0.571 splits at 1/T — 15 (compare Fig. [_J. 
The more weakly interacting particles are found predom- 
inantly in the isotropic phase, while the particles with 
larger u are partititioned across the two nematic phases. 
This tendency is observed more generally, as the case of 
two-phase I-N coexistence in the inset of Fig. [5] shows. 
Physically, it is explained easily from the form of the 
Maier-Saupe free energy: the system can lower its excess 
free energy (only) by ordering nematically. This effect 
is stronger for the more strongly interacting particles, 
which are therefore found preferentially in the nematic 
phases. 

As a consequence of fractionation, the overall parti- 
cle number fractions x a found in the various coexisting 
phases also do not vary linearly across coexistence region. 
Fig. [_] illustrates this, for a Schulz parent with the same 
polydispersity S — 0.571 as in the two examples in Fig.__] 
One notes in particular the non-monotonic variation of 
the fraction of particles in the isotropic phase before the 
second nematic phase appears. This is due to the prox- 
imity to the re-entrant I-N phase boundary in the phase 
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FIG. 4: Volume fractions x a of the various coexisting phases 
plotted as a function of inverse temperature, for a Schulz par- 
ent with polydispersity 8 — 0.571. 

diagram. 

IV. DISCUSSION 

We have studied the effect of polydispersity in ther- 
motropic liquid crystals, using a Maier-Saupe theory 
with factorized interactions. For a spread of interaction 
strengths u of around S = 0.50 of the mean and higher, 
polydispersity causes new phase splits to appear in the 
phase diagram. In particular, the theory predicts co- 
existence of two or more nematics with or without an 
additional isotropic phase. Together with this, a strong 
broadening of the I-N coexistence region is observed, to 
the extent that the phase boundary between the N and 
I-N regions becomes re-entrant within a small range of 
6. Our calculation only considered isotropic and nematic 
phases and so at least some of phase behaviour we found 
at lower T will be only metastable, the stable phases be- 
ing e.g. smectic or crystalline. Nevertheless, one would 
expect to be able to observe at least the strong widening 
of the I-N coexistence region in appropriate experimental 
systems. 

From a broader theoretical point of view, our results 
demonstrate that polydispersity in the attractive (and 
orientation-dependent) interaction between liquid crys- 
tal particles can be sufficient to cause demixing into two 
or more nematic phases. Polydispersity in short-range re- 
pulsions as caused e.g. by a distribution of particle shapes 
is not necessary, although from the results for binary 
mixtures in the Onsager model one would expect it to 
reinforce the tendency towards nematic demixing. This 



issue would merit further study but promises to be very 
challenging. Indeed, even in the polydisperse Onsager 
model without added attractive interactions it remains 
an open problem to establish the existence of nematic- 
nematic demixing at high densities. 

One may ask how crucial the assumption of factor- 
ized interactions u(l,V) = u x / 2 {l, /)u 1 / 2 (r, V) is for the 
phase behaviour that we have found. The situation is 
easiest to understand at low temperatures, where one 
has full nematic order (S(l) » 1) and the ideal (en- 
tropic) part of the free energy can be neglected. The 
system then has to minimize the excess free energy 
— \ J dl dV n(l)n(l')u(l, I'). To understand whether phase 
separation will occur, take first the simple case of a binary 
mixture, with a fraction x of particles of length (polariz- 
ability, etc) I and fraction (1 — x) of particles of length /'. 
The demixed state of two phases, each containing only 
one type of particles, has a lower free energy than the 
mixed phase if 

- l[xu(l, I) + (1 - x)u(l', I')] < -\[x 2 u(l, I) (18) 
+ (1 - xfu{l', I') + 2x(l - x)u(l, I')] 

or equivalently 

u(Z,Z') < \[u(l,l)+u(l',l')} (19) 

One can easily show that this condition in fact guarantees 
that any mixture of a finite number of components will 
separate into phases each containing only one component. 
In a fully polydisperse system, one then effects a cascade 
of demixing transitions with an ever-increasing number 
of phases as T is lowered, as is found in e.g. models of 
polydisperse copolymer blends |43| . Our observation of 
demixing into multiple nematic phases should therefore 
be generic as long as <|19H holds, i.e. as long as the at- 
tractive interaction strength between different particle 
species is always lower than the mean of the intra-species 
interaction strengths. 

In future work, one could contemplate extending the 
present study by accounting for the effects of density vari- 
ations. To a first approximation this can be done by 
dividing the interaction strengths u (l, I ') by the average 
particle volume in the given phase |20l |. In the polydis- 
perse setting, this would give an extra factor of 1/v in 
the excess free energy, where v = Jdln(l)v(l) and v(l) is 
the particle volume as a function of I. One then still has 
a truncatable free energy, but now with an excess part 
depending on two moments, S and v. Also of interest 
would be an extension to liquid crystal polymers, where 
perturbative results for small polydispersity 8 already ex- 
ist 40] and there is experimental evidence for N-N phase 
separation, at least in binary mixtures |48j |. 
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